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Abstract. In this paper we provide the characterization of all finite-dimen- 
sional Heath-Jarrow-Morton models that admit arbitrary initial yield curves. 
It is well known that afBne term structure models with time-dependent coef- 
ficients (such as the Hull— White extension of the Vasicek short rate model) 
perfectly fit any initial term structure. We find that such affine models are in 
fact the only finite-factor term structure models with this property. We also 
show that there is usually an invariant singular set of initial yield curves where 
the affine term structure model becomes time-homogeneous. We also argue 
that other than functional dependent volatility structures - such as local state 
dependent volatility structures - cannot lead to finite-dimensional realizations. 
Finally, our geometric point of view is illustrated by several examples. 



1. Introduction 

In this paper we provide the characterization of all finite-dimensional Heath- 
Jarrow-Morton (HJM) models that admit arbitrary initial yield curves. This is an 
extension and completion of a series of results obtained by Bjork et al. ji], |j| and 
H H [l0| . It is well known that affine term structure models with time-dependent 
coefficients (such as the Hull- White extension of the Vasicek short rate model jl3) ) 
perfectly fit any initial term structure. We find that such affine models are in fact 
the only finite-factor term structure models with this property, under some weak 
assumptions on the volatility structure. We also show that there is usually an 
invariant singular set of initial yield curves where the affine term structure model 
becomes time-homogeneous. This is again well known for the classical Vasicek [|6j 
and Cox-Ingersoll-Ross (CIR) Q short rate models, where the set of consistent 
inital curves is given explicitely by the model parameters. 

Practitioners and academics alike have a vital interest in finite-factor term struc- 
ture models, and the distinction of time-homogenous and inhomogeneous ones. Ac- 
cording to |l2| ] there are two groups of practitioners in the fixed income market. 

Fund managers trade on the yield curve (buy and sell swaps at different matu- 
rities), trying to make money out of it. They do not believe that all the interest 
rate market quotes are "correct" . Instead, they in general use a time-homogeneous 
two- or three-factor model, estimate the model parameters from long time series 
data, and then update the state variables (factors) each day to fit the current term 
structure. Hence the term structure is considered as a derivative based on more 
fundamental state variables (factors), such as in an equilibrium model. The dis- 
crepancies between the fitted term structure and the market prices are preceived 
as potential trading opportunities. For example, if the fitted curve is above the 
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two year and ten year swap rates, but is below the five year swap rate. Then one 
does a butterfly trade: receiving the five year rate (as one thinks it is high) and 
delivering the two year and ten year rates (as one thinks they are low compared to 
the five year rate). After this trade, one usually needs to wait for six months or 
longer for the rates to "reverse" (as predicted by the model) so that one can make 
money. Since this is a long term game, the model parameters must not change 
every day. Parameters have to be constant. If a parameter is time-varying, it is a 
factor and one needs to specify its dynamics so that one can make corresponding 
adjustments for the hedging. A state variable (factor) is time-varying, but since 
one has a stochastic model for its evolution, one can check on a daily basis whether 
its realized value lies within a statistical confidence interval or not. 

Interest rate option traders, on the other hand, often take the quoted yield curve 
data, with minimal or no smoothing, as model input. To fit the observed yield curve 
perfectly, they allow some of the model parameters to be time-inhomogeneous. 
They intend to hedge away instantly all the risks on the yield curve and only 
worry about the risk in the implied volatility structure. Yet, low-dimensionality of 
the model is desirable, since the number of factors usually equals the number of 
instruments one needs to hedge in the model. And the daily adjustment of a large 
number of instruments becomes infeasible in practice due to transaction costs. Of 
course, the model factors have to represent tradable values. But this can usually 
be achieved by a coordinate transformation. 

An HJM model for the forward curve, x i— ► r t (x), is determined by the volatility 
structure, r t <—>■ (x <—>■ a(r t ,x)), and the market price of risk. Here r t {x) denotes 
the forward rate at time t for date t + x (this is the Musiela fl3| parameterization). 
That is, the price at time t of a zero-coupon bond maturing at date T > t is given 

by 

P(t,T) = e~~ ^o r ~* rt( - x ~> dx . 

It is shown in || that essentially every HJM model can be realized as a stochastic 
equation 

dr t = ( 4~ r t + ®HJM(rt)) dt + a(r t ) dW t . _ 

\dx J (1.1) 

ro = r , 

in a Hilbert space H w of forward curves. We shall recapture the precise setup below 
in Section |[ The solution process (r t ) in general cannot be realized by a finite- 
dimensional state process. An HJM model is said to admit a finite- dimensional 
realization (FDR) at the initial forward curve r* if, roughly speaking, there exists 
an m-dimensional diffusion state process Z (factors), for some m G N, and a map 
<j) : W n — » H such that r t — <j>{Z t ). Notice that m, Z and <f> depend on r* . In 
fact, we shall be interested in those HJM models that admit an FDR of the same 
dimension at every initial curve, and this dimension is minimal in some sense. 

One of the most basic examples for a finite-dimensional HJM model that fits 
any initial curve is the Hull- White extended Vasicek model |l3| for the short rate 
R t :=n(0), 



dRt = (b(t) - (3R t ) dt + pdW t , Ro = r*(0). 
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Here W is a one-dimensional Wiener process defined on some stochastic basis 
(f2, (J-t), J 7 , Q), where Q is the risk-neutral measure. This example will be com- 



pletely recaptured in Section 6.1, where also all the subsequent functions are ex- 
plicitly given. The coefficients (3 > and p > are constant, and b(t) is a function 
determined by the initial forward curve r* . The corresponding volatility structure 
is constant, 

a{r,x) = u{x)= pe- px , (1.2) 

and the forward curve has an affine dependence on R t , 

r t (x) = A H wv(t, x) +B HWV (x)R t . (1.3) 

Whence a one-factor affine term structure model with time-dependent coefficients 
(strictly speaking, the realization is given by the two-dimensional process (t,R t )). 
For those inital curves for which b(t) = b we obtain a time-homogeneous process 
(Rt) and Ajjwv(t, x) = Ahwv{ x )- in this example we also can easily perform a 
coordinate transformation that leads to a different state process Z. Indeed (as in 
d Section 4.1]) let 

dZ t = -/3Z t dt + p dW t , Z = 0. 



Then we have 



and hence 



R t = Z t + e~ f3t r* (0) + / e^'^^s) dt, 



r t (x) = A HWV (t,x) + B HWV (x)Z t , (1.4) 

for some function Ahwv ■ Thus we obtaine d a simpler state-process Z with a 

similar (affine) functional form of r t (x) as in (|l.3| ). But Z is not a tradable value 

and hence cannot be directly used for hedging. 

There is a substantial literature providing sufficient conditions for the existence 

of finite-dimensional HJM models (see e.g. || for further reference). A systematic 

study from a geometric point of view has been made by Bjork et al. (l[ [|, §], and 

§ | 0, see also @. In § Bjork and Svensson give necessary and sufficient 

conditions for the existence of FDRs. Their key argument is the classical Frobenius 

theorem. Therefore they define a Hilbert space, 7i, on which d/dx is a bounded 

linear operator. This space consists solely of entire analytic functions. It is well 

known however that the forward curves implied by a CIR short rate model are of 

the form r t (x) = A C ir{x) + Bcm{x)R t , where 

e ax _ ^ be ax 
A C ir(x) = d and B CIR {x) = (1.5) 

e + c (e ax + c) z 



for some a, b, c > and d > (as will be recaptured in Section p.2[ ). But Aqir and 
Bcir are not entire analytic functions since their extensions on C have a singularity 
at x = (lnc + in) /a. Hence the CIR forward curves do not belong to H, whence 
the Bjork-Svensson || setting is too narrow for the HJM framework. 

In [[l0| we recently succeeded to overcome this difficulty. The aim of the present 
article is to bring these results to the attention of the financial community. We 
shall make clear why, under fairly general assumptions, finite-dimensional HJM 
models are necessarily affine, and provide a deeper understanding of some important 
examples. 
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The remainder of the paper is as follows. In Sections [2| and [5| we give the precise 
setup and definition for a (local) HJM model and an FDR around an initial curve, 
respectively. In Section || we recapture the main results (without proofs) from [[To) 
and We only consider functional dependent volatility structures, which in- 

cludes essentially all interesting examples. This is also justified by Section |5|, where 
we show that a local state dependent volatility structure cannot admit an FDR, un- 
less it is constant. In Section ^ we illustrate our results by explicit calculations for 
two-dimensional HJM models, which turn out to be either of Hull- White extended 
Vasicek or CIR type. 

2. HJM Models 

Let (f2, (Ft), T ', Q) be a filtered probability space satisfying the usual conditions, 
and W = (W 1 , . . . , W d ) a d-dimensional Wiener process, d > 1. We suppose that 
Q is the risk neutral measure for the subsequent models. Indeed, for the existence 
of an FDR it is irrelevant whether we are under the physical measure P or under 
<Q> ~P (see also §, Remark 7.1.1]). 

We follow H , where r t is regarded as an element in the Hilbert space of forward 
curves H w . This space consists of absolutely continuous functions h : R>o — > M, 
and is equipped with the norm 

\\h\\l := \h(0)\ 2 + f \d x h(x)\ 2 w(x)dx, 

JWL> 

where w : M>o — > [l,oo) is a non-decreasing C 1 -function such that w" 1 / 3 is inte- 
grable on R>o- For example, w(x) — e ax or w(x) = (1 + x) a , for a > or a > 3, 
respectively. The closed subspace H w q of H w , consisting of functions h £ H w with 
lim^^oo h(x) = 0, has the property that 

(9, h)^g J h, 

where J h denotes the definite integral x i— > h(y) dy, is a continuous bilinear 
operator from H w $ x H w q into H w . Hence whenever a = (a\, . . . , <7d) is a (bounded 
and) locally Lipschitz continuous map from H w into H!~ then the HJM drift term, 

d - 

a hjm (h)=J2^(h) Vi(h), (2.1) 
i=i 

is a well-defined (bounded and) locally Lipschitz continuous map from H w into H w . 

On H w the (unbounded) closed operator A — d/dx generates the strongly con- 
tinuous semigroup (St) of right-shifts 

S t h = h(-+t), t>0. 

The same holds for the restriction Aq of A to H w ^. It can be shown that the 
domain of A is D(A) = {h 6 H w n C^R^o) | d x h G i/^} and similarly D(A ) = 
{/i G H wfi n C^R^o) | 9 x /i G #™,o}- We consider A 2 = Ao A and all higher order 
powers of A and Ao. By induction, we have 

D(A°°) := p| D(A n ) = {he H w n C°°(M> ) | («9 x )"/i G Vn > 1}, 

n>0 

C(AS°) := f| D(A^) = {/i G £T Wi0 n C°°(R> ) | (d x ) n h £ H w>0 , Vn > 1}. 

n>0 
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It is clear that D(Ag°) c D(A°°). Equipped with the sequence of seminorms 
Pn{h) = XT=o II^^ILi n > 0, D(A co ) and D(A§°) become Frechet spaces, and 
A acts as a bounded linear operator on them. We refer to Q for the theoretical 
background. In fact, D(A°°) and D(A^) are in some sense the largest subspaces 
of H w and H Wt o, respectively, with this property. The functions from (fL3]), (1.4) 
and (|J), A HWV (t, •), A HWV (t, •), A C ir, B hwv , B C ir, all lie in D(A°°), the last 
two even in £>(^4g°). 

By a solution (r t ) of (1.1) we shall always mean a continuous mild solution (see 
(5[ ||). That is, an H w -vahied continuous adapted process (rt) which satisfies 

r t = S t r*+ / S t - s a H jM(rs)ds + y] / 5 t _ s a,(r s ) <W S \ (2.2) 

In the classical HJM notation, where f(t, T) = r t (T — t), this becomes the familiar 
expression 

f{t,T) = f{0,T)+ / a H jM(8,T)da + J2 / ^{s,T)dW:, < t < T, 
Jo i=1 

where a HJM {s,T) := a HJM (r s ,T - s) and 6\(s,T) := cTi(r s ,T - s). If ( p^ ) only 
holds locally for i replaced by t A r, for some stopping time r > 0, then (r t ) is a 
Zoca/ solution. 

Since the volatility structure cr = (eri , . . . , aj) determines cihjm by (2.1), the 
following terminology is justified. 

Definition 2.1. Let U be a convex open set in H w . A (local) HJM model in U is a 
map a = (o"i, . . . , ad) ■ U — > swc/i that (1.1) admits a unique U-valued (local) 
solution for every initial curve r* G U. 

It is shown in |)], Section 5.2] that a :U — > ffj* is a local HJM model in U if it 
is locally Lipschitz continuous. The restriction to U is convenient since it allows to 
incorporate important examples such as the CIR model, where U is the half space 
{h€H w \ h(0) > 0}. 

3. Finite-Dimensional Realizations 

In this section we give the rigorous definition of an FDR as it was sketched in 
the introduction. 

Let U be a convex open set in H w and a a local HJM model in U. Let Tq G 
UC\D(A ca ) and neN. 

Definition 3.1. We say that a admits an n-dimensional realization around r$ 
if there exists an open neighborhood V of 7q in U n D(A°°), an open set U in 
R> x W 1 - 1 , and a C 00 -map a : U x V -> U n 13 (A 00 ) swc/i tfwrf 

i) r £ a(J7, r) /or a/Z r 6 

ii) D z a(z,r) : R n — ► D(^4°°) is injective for every (z,r) G J7 x V , 

iii) a(zi,ri) = a(z 2 ,r 2 ) implies D z a(z\, ri)(R n ) = D z a(z 2 , r 2 )(R") /or (zj, r,) g| 
UxV, 

iv) /or every r* £ V t/iere exists a U-valued diffusion process Z and a stopping 
time t > such that 

rtAr = a(Z tA T,r*) (3.1) 



is the (unique) local solution of (1.1) with ro = r* 
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Thus we only consider FDRs in U n D(A°°). This seems to be a restriction 
first, since the original local HJM model is defined on U. However, we shall see in 



Proposition 4.3 below that in most interesting cases the FDRs are necessarily in 
U n D{A°°). Also we should not worry about U C M> x It will become 

clear that the first component of the diffusion Z can always be chosen to be the 
time t > 0. 



Definition 3.1 states that an n-dimensional realization around 7q implies the 
existence of an FDR at every point r* in a neighborhood of Tq , and these FDRs have 
a smooth dependence on r* . In fact, by i) and ii), each a(-, r*) : U — > U n D(A°°) 
is (after a localization) the parametrization of an n-dimensional submanifold with 
boundary, say M. r * , of U H D(A°°), and (|3.1|) says that 



r tAT e Mr', Vi > 0. 

Condition iii) implies that two such leafs A4 ri and M r9 can only intersect at points 
where their tangent spaces coincide. According to |10|] , the family {AlVlrev is 
called an n-dimensional weak foliation on V. 

The existence of FDRs around a point is assured by an extended version of the 
Frobenius theorem (JToj] ) on the Frechet space D(A°°). The Frobenius theorem has 
also been used by Bjork et al. |5|, |) on the Hilbert space H, which however has the 
drawbacks mentioned in the introduction. 

A striking feature of the Frobenius theorem is that it brings together an algebraic 
condition (dimension of a Lie-algebra) with the analytic problem of the existence 
of an FDR (weak foliation). The former condition can in many cases be explicitly 
checked, as it is exemplarily carried out in Section 0. We do not intend to go further 
into the theory of Lie-algebras here, but refer the interested reader to |l(J. We only 
mention that, under the appropriate assumptions (in particular, a : U (~l D(A°°) — > 
D(A^) d has to be C°°), a necessary condition for the existence of an n-dimensional 
realization around is that the Lie-algebra generated by the vector fields 

1 d 

H(r) := Ar + a HJM (r) - -22 D(X i( r ) a i( r ) ( 3 - 2 ) 

i=i 

and <7i, . . . , ad onWn D{A°°) is at most n-dimensional. This algebraic condition 
typically yields an obstruction for the examined model to admit FDRs, as it is 
shown in Section |^. 



4. Functional Dependent Volatility 

In this section we have the idea of a being sensitive with respect to linear func- 
tionals of the forward curve. That is, we let a : U — > H d be a local HJM model 
in U and suppose in addition that tJi(r) = <fri(£i(r), . . . ,£ p (r)), for some p G N, 
where 4>i '■ O C MP — > D(A^) are smooth maps and £%, . . . , £ p are continuous linear 
functionals on H w . This includes essentially all interesting examples. Indeed, we 
may think of 

1 f xt 

£i(r) = — / r(y) dy (benchmark yields) 

x i JO 

or 

li(r) — r(xi) (benchmark forward rates), 
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for some fixed tenor x\ , . . . , x p > 0. For a short rate model we simply have p — 1 and 



£(r) = r(0) (see (1.2) for the trivial constant case £(r) = 1 and (f>(z) = pexp(— /?■ )). 
This idea is formalized by the following regularity and non-degeneracy assumptions 

(Al): <Ti = fa o £ where I G L(H W ,W), for some p G N, fa : O -» £>(A§°) are 

C°°-maps, for 1 < i < d, and is an open set in R p containing £(U). 
(A2): 4>i(z), . . . , (j>d(z) are linearly independent, for all z G O. 
(A3): The linear map (£,£oA,...,£oA q ) :L»(A°°) -> Rp(?+i) 

is open, for every 

finite q > 0. 

Assumptions (Al) and (A2) are clear, only (A3) needs some further explanation. 
Intuitively, it says that the following interpolation problem is well-posed on D(A a °): 
given a smooth curve g : R>o — -> M, for any finite number of data of the form 
1 = (£(g),£(d x g),...,e((d x ) q g)) G we can find an interpolating function 

h G D(A°°) with (^(/i), ...,£ o A q (h)) = X. Hence degenerate examples such as 
p = 3 and £{h) = (/i(0), /i(l), J /i(x) dx) are excluded. Indeed, here we have 
£ o A(h) = (d x h(0), d x h(l), h(l) - h(0)). Thus the rank of (£, to A): D(A°°) -» M 6 
is 5, whence this map is not open. Assumption (A3) is always satisfied for short 
rate models, since there p = 1 and £(h) = h(0). 

Example 4.1. Constant direction volatility. Let Ai, . . . , be fixed linearly inde- 
pendent vectors in D(A™), and 

d 

<7j = 4 o £ with <j>i{z) = 2 G 0, (4.1) 

i=i 

for some C°°-map $ = ($y) on O with values in the regular d x d- matrices. This 
example has been extensively studied by Bjork et al. || |). 



The vital importance of this example will be made clear by Proposition 4.4 and 
Theorems 4.6 and 5.1 below. 

Now let a satisfy (A1)-(A3). This implies in particular that a : U — > is 
locally Lipschitz continuous and hence a local HJM model in U. 

Let rj gMfl D(A°°) and neff. There is a lower bound for the dimension of an 
FDR around rjj. 

Proposition 4.2. Suppose a admits an n-dimensional realization around Tq. TTien 
n > d+ 1. 



Proo/. This is flO, Proposition 4.8]. □ 



Hence there exist no one-dimensional realizations around 7g, even if d = 1. This 
confirms the well-known fact that a short rate model (R t ) that fits every initial curve 
r* in a neighborhood of 7q in Z)(A°°) contains necessarily some time-dependent 
parameters. Such that Z t — (t,R t ) yields a 2-dimensional realization around rg. 
This also proves a conjecture in || (see Remark 7.1 therein). 

Since low-dimensionality of the state process is preferred, it would not make much 
sense to look for n-dimensional realizations with n > d + 1 (although this can be 
done under some non-degeneracy assumptions, as it is carried out in fTc[] ). Rather 
we restrict our attention to minimal (that is, (d -I- l)-dimensional) realizations 
around Tq. In this case our focus on FDRs in D(A°°) turns out to be no restriction 
at all. 



8 



DAMIR FILIPOVIC AND JOSEF TEICHMANN 



Proposition 4.3. Suppose M. is a (d+1) -dimensional C°° -submanifold (with bound-^ 
ary) o}U that is locally invariant for (|l.l[). That is, for every r* G Ai, the U-valued 



local solution (rt) of (1.1) satisfies 

r tAT G M, Vt > 0, 

for some stopping time r > 0. Then necessarily Ai C D(A°°), and Ai is a C°° - 
submanifold (with boundary) of D(A°°). 

Proof. This is (TiJ Theorem 3.2]. □ 

We now can cite jn| Theorem 4.10]. 

Proposition 4.4. Suppose a admits a (d + I) -dimensional realization around Tq. 
Then a is of the form d4.1|) (constant direction volatility) on V, where V is given 



by Definition 3.1 



So far our considerations were local. Now let U denote the set of all 7q G 
U (~l D(A°°) around which a admits a (d + l)-dimensional realization. Is it pos- 
sible that U = U n D{A°°)1 In general the answer is no. Indeed, suppose U is 
connected. Then, by a simple continuity argument, a is of the form ( [hip with the 
same Ax, . .., A<j everywhere on U. It then follows from |l^, Theorem 4.10] that U 
must not intersect with the singular set 

S:= {fteWn D(A°°) I v(h) G (Ax, . . . , X d )}, (4.2) 

where 

v(h):=Ah + a HJ M(h) (4.3) 

and (Ax, ■ ■ ■ , Ad) denotes the d-dimensional subspace spanned by Ax, • • ■ , A<j. The 
set E is not empty in general, but it is small. 

Lemma 4.5. The set S is closed and lies in a subspace G of D(A oc ) with (d+1) < 
dimG <d 2 +d+l. 

Proof. This follows from jl0| Lemmas 4.5 and 4.11]. □ 

We now can summarize and state the main result, which says that only afhnc 
term structure models admit a minimal FDR around any initial curve, see [ flO[ 
Theorem 4.13]. 

Theorem 4.6. Suppose there exists an open connected set U in U n D{A oc ) such 
that a admits a (d + 1) -dimensional realization around every r^ G U. Then there 
exist linearly independent vectors Ax, • • • , Xd in D(A^) such that a is of the form 



(4.1) on U , and U fl S = 0, where E is given by (|L2j). Moreover, for every r* G U 
there exists e > and a C 00 -map *(-,r*) :[0,t)^Wn D(A°°) with *(0, r*) = r* , 
an M d -valued diffusion process Z with Zq = and a stopping time < r < e suc/i 

r tAr = vI/(iAT,r*) + ^Z t l Ar A l (4.4) 

i=I 

is t/ie U-valued local solution of (|l . l| ) . Moreover, ^(-,r*) is £/ie unique solution to 
the evolution equation in D(A°°), see ( |3.2| ), 

^u(t) = /*(«(*)), u(0)=r*. (4.5) 
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//, in addition, U = U D D(A°°) \ X then for every r* 6 X i/iere exists an 
Mr -valued time-homogeneous diffusion process Z with Zq = cmd a stopping time 
< r < e smc/i i/iai 

d 

r, Ar =r*(tAr)+^Z| AT A i (4.6) 



invariant for (1.1). 



is the U -valued local solution of (1.1), and r tAT G S /or allt > 0. Hence £ is locally 



Thus here we have the announced decomposition of the space of initial for- 
ward curves U n _D(A°°) into t/, where the (affine) factor model becomes time- 
inhomogeneous, and the singular set S, which is invariant for the model dynamics 
and where the (affine) factor model becomes time-homogeneous. This phenomenon 
was known for some particular models, and it now is proved in full generality. We 
shall further illustrate this result for the case d = 1 in Section IBl 



Remark 4.7. The vector field /i in (4.5) can be replaced by v and the statement 



about the affine form (4.4) remains true if r is chosen such that r tAr is U HU- 



valued. This follows since fJ-(h) — v{h) G (A l7 . . . , A^) for all h € U Dli (= U if 



U = U\ Yi), which yields a straightforward modification o/\I/ and Z in (4.4) 



5. Local state dependent volatilities 

Continuing the discussion at the end of Section |[ we are formally given smooth 
vector fields Xx, —,X m onMfl D(A°°). Involutivity - the algebraic condition of 
the Frobenius Theorem - amounts to saying that for any i,j = 1, ...,m there exist 
smooth functions Ai, A m : Uq — > R, such that 

m 

[Xi,Xj]ih) = ^Xi(h)Xi(h) 
i=i 

for h in some smaller domain Uo C U n D(A°°). We recapture some facts, which 
have been derived in jio) . The calculation of Lie brackets can be performed on the 
space C°°(]R>o,K), since D{A°°) is a subspace of C°°(R>o,K) and the vector fields 
that we consider are defined on the latter. Lie brackets are calculated by Frechet 
derivatives 

[X t , XA{h) = DXi(h) ■ Xj(h) - DX,(h) ■ Xi(h), 

and calculating Frechet derivatives is equivalent to derivation of _D(A°°)-valued 
curves 

DXi(h)-v= j t \ t =oXi{h + tv). 

We want to demonstrate that local volatilities cannot lead to finite-dimensional 
realizations, since the algebraic conditions are not satisfied. For simplicity we shall 
assume that d = 1. Given <fi : ]R>o x M — > M, we define <r(h)(x) — 4>(x, h(x)) for 
h G U n D(A°°) and x > 0. By cartesian closedness (see Q) we obtain for the 
Frechet derivative 

(Do-(h) ■ v)(x) = ^-\ t=0( j)( x , h(x) + tv{x)) 
at 

= 4>'{x, h{x))v{x) 
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for any x > 0. Consequently D<r(h) ■ v = ((f)' o h)v, where we have an ordinary 
multiplication of functions. With o we denote the composition with respect to the 
second variable and with .' we denote derivation with respect to the second variable. 
Therefore 

Kh) = + (<f> o h)( J <t> o h) - \{{<t>y o h)(4> o h), 

see (ft.2|). Conditions on (f> can be found such that all parts of the drift vector field 
and the local volatility are smooth maps. 

Theorem 5.1. Assume that [p, o~](h) 6 ((j(h), jj,(h)) for h in some open set Uq C 
IA (~l D(A°°). Then o~(h)(x) = pexp(—f3x) for h £ Uq, for some constants p and (3, 
which is the Vasicek volatility structure. 

Proof. We have to calculate one Lie bracket: 

D^(h)-v = ^-v + (U' oh)v)(( 4> o h) + U o h)( ( (6'oh)v)- 
dx J J 

-~(0" oh)((f)oh)u -htf ohfv 

for h € UC\D(A ca ) and v £ D(A°°). The derivative with respect to the first variable 
of 4> is denoted by d\ . This leads to 

[fJ,,a](h) = {<p' oh)-^h+(d 1( f> oh) + (<t>' oh)(<p oh) (^J (poh\ + 

+(<poh) (^J ((j)' oh){(j>oh) \ - ^((f)" O h)((f) O h)((f) O h) - ^((f)' O h) 2 ((f) O h) - 

o h)^-h - (cf)' o h)((f) o h) (^J ((f) o h)j + ^((f)' o h)((f)' o h)(<f> o h) 
= (di<j) oh) + ((f)oh) (^J ((f)' o h)((f) o h)j - ^((f)" o h)((f) o hf. 

We now shall evaluate the equation 

\H,<r](h) = \(h)a(h) + \ Q (h)p(h) 

for ft e Mo. Assume that Ao(ftoo) ^ for some h^ G Uq. We can isolate 4-h, then 
we take a sequence (h n ) n >o m Uo with h n — > ftoo in D(A m ) as n — > oo but not 
in £>(A m+1 ), where m is chosen such that \o(h n ) — > Ao(/ioo) (this is possible by 
applying Theorem 2.3 of and adjusting Uq). Hence the only term which does 
not converge up to order m is -^h n . Consequently Ao = on Uq. Next we analyze 
the resulting equation 

[fi,a](h)=X(h)a(h), 

where we proceed by the same reasoning. 

First we assume that we are given a point ho G Uo and xo > such that 
(f)(xo, ho(xo)) ^ 0, then we can divide by a(h)(x) for h in a neighborhood of ho and 
i in a neighborhood of xo and identify a logarithmic derivative. To this equation 
we apply again the operator and obtain 

(Of In \(f)\ oh) + (<9i(ln |0|)' o h)h' + ((f)' o h)(4> oh)- -(di(f)" o h)((f> o h) - 
-~(4>'" o h)h'((f) oh)- -(</>" o h)(di(f) oh)- -((f)" o h)((f>' o h)h' = 0. 
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We can again isolate h! = 4-h, which leads to a contradiction as before, therefore 
its coefficient has to vanish identically. This leads to the following two equations: 

(0? In |0| oh) + (0o/i)(0' oh) - -(9i0" o/j)(0o/i) - (0" oh){di<t>oh) = 

and 

(0i (In |0|)' o h) - ~(0"' o /i)(0 o h) 2 - (0" o fi,)(0' o h) = 0. 

We take these two equations and evaluate them for h in a neighborhood of ho and 
a; in a neighborhood of xq, which leads to the equations 

0i 2 ln|0|+00'-i0 1 2 (0"0)=O 

0i(m|0|)'-l(0"0)' = O. 

Taking the derivative .' in the first and 0i in the second and finally the difference 
of the resulting equations we obtain 

(00')' - 

and therefore 4>(x,y)-j^4>{x,y) — g(x) for (x,y) in a neighborhood of (xq, ho(xo)) 
with some smooth function g : M>o — > K. This equation has a smooth solution if and 
only if J^(f)(x, y) = 0. Therefore 2 In |0| = 0, which leads to <f>{x, y) = pexp(—(3x) 
for (x,y) in a neighborhood of (xq, ho(xo)). By continuity we obtain the global 
result. □ 



Remark 5.2. The same method can be applied to the d-dimensional case. This 
leads to a similar assertion, however, the calculus is rather ambitious. 



6. Characteristic Examples 

In the seminal papers and finite-dimensional realizations, in particular the 
Hull- White extensions of the Vasicek and CIR-model, are considered for the first 
time from the geometric point of view. In addition to their excellent treatment 
(compare Section 5 of Q or Section 7 of 0), we prove that the Hull- White exten- 
sions of the Vasicek and CIR model are the only 2-dimensional local HJM models 
and we demonstrate the importance of the corresponding singular sets. The same 
type of analysis can also be performed in higher dimensional cases, which will be 
done elsewhere. At the end of this section we provide an example of how to em- 
bed the Svensson family as a leaf of a weak foliation associated to a functional 
dependent volatility structure. 

We let again d = 1. Starting with a functional dependent volatility structure a 
and U as in Section Wand assuming the existence of 2-dimcnsional realizations on 



U n D(A°°) \ E (see ( |4.2| )), we necessarily arrive by Proposition 4.4 at a constant 
direction volatility on U. We shall show that this volatility is either of the Vasicek 
or CIR type. 

In view of (A2) we have a ^ on U, hence we can write <r(r) = 4>{r)\ for 
r e U, for some A € D{A^) \ {0} and a smooth map : U — > E, such that without 
loss of generality > (by a slight abuse of notation, the meaning of here is 
different from Section ^) . We want to specify under which conditions this volatility 
structure admits 2-dimensional realizations and how they look like. This is already 
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done in Section 7.3 of ||, however, their special setting does not allow to treat the 
CIR-case. 

Writing ip(r) := 4>{r)(D(f>{r) ■ A), we obtain for r eWn D(A°°) 
Da(r) ■ h = (Dcj)(r) ■ h)X 
Da{r) ■ cr(r) = cf>{r){D^{r) ■ X)X = ^(r)X 

H{r) = ^r + 4>{r) 2 \ J A - ~V(r)A 

£>/z(r) • h = —h + 2<t>{r)(D<f>(r) ■ h)X / A — -(DibM ■ h)X. 
ax J 2 

Consequently we can calculate the Lie bracket 

[ M , <r](r) = 0(r)^A + 20(r)V(r)A J X- l#r)(Ity(r) ■ A)A - 

d (' 1 

- (U0(r) • — r)A - 0(r) 2 (Dcb(r) ■ X J X)X + -V(r)(X>^(r) • A) A. 

We assume [/z, cr](r) € (A) on Ur)D(A oa ), which follows from the Frobenius condition 
and is justified by Lemmas 2.12 and 3.4 of [|l0). We can divide by cj>{r) and obtain 
an equation 

— X + 2ib(r)X [ X-6(r)X = 
ax J 

with some smooth function 9 : U n D(A°°) — > R. There are consequently two cases: 

i) If A and A J A are linearly independent in D(A°°), then by derivation with 
respect to r we obtain that tp and are constant, say 2ip(r) = a and 0{r) = b 
with real numbers a and b. Defining A := J X we obtain finally a Riccati 
equation for A, which yields the CIR-type if a ^ or the Vasicek-type if 
a = 0: 

-^-A + ^A 2 + 6A = A(0), A(0) = 0. (6.1) 
ax 2 

The Ho-Lee model is considered as particular case of the Vasicek model for 
6 = 0. 

ii) If A and X J X are linearly dependent in D(A°°), then we necessarily obtain an 
equation of the type 

— A + 6A = 0, 

ax 

which yields that A is vanishes identically, since otherwise A and X J X are 
linearly independent. This case was excluded at the beginning. 

Notice that by ( |6.1[ ), A(0) = if and only if A = 0, which is not possible. Hence 
a fortiori we have A(0) ^ 0, such that by rescaling we always can assume that 
A(0) = 1. This observation slightly improves the discussion in Section 7.3 in 0. 

By the definition of ip we have Dcf> 2 (r) ■ X = a, hence we obtain the following 
representation for 0. We split D(A°°) into RA + E, where E := ker evo. We denote 
by pr : D(A°°) — > E the corresponding projection. Then 



4>{h) = ^aev (h)+r,(pr{h)), (6.2) 

where r\ : pr(U OD(A CO )) C E — > R is a smooth function (compare with Proposition 
7.3 of fl). 
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Recalling (4.2), we have 

heUn D{A°°) | v{h) = Ah + 4>{h) 2 X I A £ (A) 



Thus, if A and A J A are linearly independent in D(A°°) then any ft, £ £ is necessarily 
of the form 

ft = ai + a 2 A 2 + a 3 A 

in all cases for some real numbers a^. By the particular representation of <j> we 
obtain that a% = aa% + 3(03), where g is some smooth real function derived from 

aa\ + ■q{a 2 h 2 + a 3 A) = a 2 . 

By Fl we denote the local (semi-)flow of a vector field X on U (~l D(A°°). The 
leaves through r* of the weak foliation are given by the local parametrization 

d 

dx 



(u ,u 1 )^Fl^ o (r*) + u 1 —A 



if r* does not lie in the singular set £. If r* £ £, then the leaf is a one dimensional 
immersed submanifold of (1,A,A 2 ). Notice that the stochastic evolution of the 
factor process takes place in the ui-component, see Theorem 4.6 and Remark 4.7. 
We summarize the preceding results in the following theorem. 

Theorem 6.1. Let a and IA be as above. Assume that a admits a 2-dimensional 
realization around any initial curve r* £ hi H D{A°°) \ S. Then there exists A £ 
D^Aq 3 ) and a function cj) :IA —* M>o of the form ( |6.2| ) such that o~(h) = <j>(h)\. The 
singular set S is a (possibly empty) subset of (1, A, A 2 ), where A = J A satisfies the 
Riccati equation (6.1). The local HJM model is an affine short rate model. That 
is, for every initial curve r* £ VI PI D(A°°) there exist functions b : R>o x R — > M, 
9 : R>o — > R, ^4 : [0, e) x R>o — > R anrf a stopping time r > smc/i i/iat 

r tA rW=#AT,j;) + A(j:)J? iAr (6.3) 



is f/ie unique U -valued local solution to ( |l.l| ) anrf i/ie s/iori raies i?t = r t (0) follow, 
locally for t < t, a time-inhomogeneous diffusion process 



dR t = b(t, R t )dt+ s/aRt + 9{t) dW t . 

This process becomes time-homogeneous if and only if r* £ S, and then r tAT £ S 
for all t > 0. 



Proof. We know that A(0) / 0. Hence follows from (|4.4|) . The rest of the 

theorem is a consequence of Theorem |4.6| and the preceding discussion. □ 

6.1. The Hull- White extension of the Vasicek model. We consider the vola- 
tility structure of the Vasicek model: a(r){x) — pexp(— (3x) = p\ with p > and 
f3 > 0, for r £ U n D{A°°) = D(A°°) and x > 0. Then by the above formulas 

[Mi ct] = -0pX. 
The singular set E is characterized by 

4-ft + ^77 exp(— — exp(— fix)) — cexp(—/3x) 
dx p 
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for some real c. Therefore <Z2 is some fixed value, namely 02 = and 01,03 are 
arbitrary. Consequently the singular E set is an affine subspace for the fixed values 
p,(3: 

P 2 

h = ai - y A 2 + a 3 A. 
Going back to traditional notations for the Vasicek model we write 

A(x) = I(l-exp(-/3x)) 
B v {x)=A'(x) =e~ Px 
A v (x)=bA(x)-^-A(x) 2 , 

then h lies in the singular set E if and only if 

he A v + (B v ) 

for some value b (which becomes an additional parameter in the short rate equation) . 
The solution for r* in the singular set reads as follows 

r t =A V + B v R t 
dR t = [b-0Rt)dt + pdW t , 

where Rt = evo(r t ) denotes the short rate, which is the Vasicek short rate model. 

Outside the singular set E we have a 2-dimensional realization. First we calculate 
the deterministic part of the dynamics 

r° P 2 

Fl u ( r *)(x) = S Uo r*(x) + J S Uo - s (—exp(-0x)(l-eDq)(-l3x))ds 

o 2 f u ° d 

= S uo r*{x) + y J d^ A ) 2(x + u o- s ) ds 

2 2 

= r*{x + w ) + yA(a; + u ) 2 - yA(i) 2 . 

If we identify uq with the time variable t, which is possible since the stochastics only 
occurs in direction of Bhwv (see Remark 4.7), we obtain by direct calculations for 
(0) 

2 2 

r t (x) =r*(x + t) + ^A(x + t) 2 - ^A(x) 2 + A'(x)Z t 
dZ t = -j3Z t dt + pdWf 

A parameter transformation yields the customary form, namely 

t 



Rt = e~ 0t r* (0) + / e-W-^bis) ds + Z t 
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This yields the following expressions: 



2 2 
P Kf- , ^2 P \ I N2 



Ahwv (t, x) = r*(x + t) + !jA(x + ty - '— A{x) 

- (A'(x)) 2 r*(0) - A'(x) / e-^- s) b{s) ds 



B HWV {x) = B v {x) = A'(x) 

dR t = (b(t) - 13 R t ) dt + pdW t 
Tt = AHWv(t) + BjiwvRt 

Ht) = j t r*(t) + f3r*{t) + ^(1 - exp(-2/3t)). 

The functions Ajjwv an d Bjjwv are solutions of time-dependent Riccati equations 
constructed by geometric methods. The equation for b follows from the fact that 
A HWV (t,0)=0. 

6.2. The Hull- White extension of the CIR model. We proceed in the same 
spirit: a(r) := p^J evo(r)X for p > 0. The volatility structure is defined on the 
convex open set U = {euo(r) > e} for some e > 0. The function A := J A satisfies 
(in certain normalization) a Riccati equation, namely 

A A + ^A 2 + /3A = 1, A(0) = 0. 
dx 2 

We obtain the solution (see e.g. H Section 7.4.1]) 



2exp(z v //32 + 2p 2 )-l 



- I3)(exp(x^/f3 2 + 2p 2 ) - 1) + 2^ + 2p 2 ' 

Under this assumption we can proceed as above: The singular set S is determined 
by the equation 

d 



hence 



u(h) = ^-h + p 2 ev (h)AA' e (A) 
dx 



P 2 

h = a 1 + ^-aiA 2 + a 3 A. 



Again a\ and 03 can be chosen freely, which completely determines S. Traditionally 
one writes the singular set in the following form: 

Acir = bA 

Bern = 1 - PA ~ yA 2 = A' 

with some additional parameter b and we obtain equally that h lies in £ if and only 
if 

h e Acir + {B C m) ■ 
The short rate dynamics follows the known pattern: 

Tt = Acir + BciRRt 
dR t = {b - (3R t ) dt + pyJ~R t dW t 
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for r* s S. Outside the singular set we have a 2-dimensional realization. First we 
calculate the deterministic part, by the variation of constants formula, 

Fi: (r*)(x) = S Uo r*(x)+p 2 / ^(r*)(0)(5 Uo _ 5 (A'A))(x) ds. 

Jo 

Identifying uq with the time parameter yields the following formula 2-dimensional 
realization, which is derived by direct calculations, 

r t = Fl»{r*)+k'Z t 

dZ t = -/3Z t dt + P y/c(t) + Z t dW u 

where c(t) = Fl%(r*)(0). The short rate is given through R t = c(t) + Z f and 

dRt = (/3c' (t) - (3Z t ) dt + py/R~ t dW t 

= (b(t) - f3R t ) dt + py/R~ t dW t . 

Notice that A(0) = A'(0) = 1 by the Riccati equation and b(t) = d(t) + (3c(t). 

This formula closes the circle with the classical Hull- White extension of the 
CIR-model: 

A H wci R (t,x) = Fl u t (r*)(x) ~ c(t)A'(x) 
Bhwcir — Bcir = A' 

dRt = (b(t) - (3R t ) dt + p^fR t dWt 

rt — AjiwciR{t) + BnwciRRt 
6 (i) + 

c(t) = r*(t) + p 2 f c(s)(AA')(t - a) ds. 
Jo 

Again this is a geometrical construction of solutions of time-dependent Riccati 
equations. 

6.3. Fitting procedures as leaves of foliations. A popular forward curve- 
fitting method is the Svensson Jlj| family 

Gg(x, z) — z 1 + z 2 e~ z -> x + z 3 xe~ ZsX + z 4 xe~ ZsX . 

It is shown in (7| that the only non-trivial interest rate model that is consistent 
with the Svensson family is of the form 

r t = Z\gx + ■■■ + Z?g 4 , (6.4) 

where 

3i (x) = 1, g 2 (x) = e' ax , g 3 (x) = xe~ ax , g 4 (x) = xe~ 2ax , 
for some fixed a > 0. Moreover, 

Z\=Zl Zf = Zle- at , Zf = Z*e- 2 « t (Z 4 > 0) 
and Z 2 satisfies 

dZ 2 = (Zf + Zf - aZf) dt + yfaz} dW t . (6.5) 
Here W is a real-valued Brownian motion. 
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We now shall find a 2-dimensional local HJM model that is of the form ( |(i.4| ) 
whenever tq = JZj—i Zj9j with Z4 > 0. In view of fl6.5p , a candidate for a is given 
on the half space U := {£ > 0} by 



a(h) = y/a£(h)g 2 , 

where £ is some continuous linear functional on C(R>o,R) with £(g\) = £{()i) — 
^(53) = an d ^(54) = 1- Straightforward calculations show, for h G U (~l D(A°°), 

v(h) = Ah + £(h)g 2 ~£(h)gi 

\M,a]{h) = ~a^d(h)g 2 - ^=^=92- 

2y/a£(h) 

(the clue is that £ o a = 0) . Hence indeed the Lie algebra generated by a and has 
dimension 2onWn Z^A 00 ) \ S. 
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